Oscillations in the inflaton potential? 

Cedric Pahud^^^, Marc Kamionkowski^, and Andrew R. Liddle^ 



OO 

o 
o 



^Astronomy Centre, University of Sussex, Brighton BNl 9QH, United Kingdom 
^California Institute of Technology, Mail Code 130-33, Pasadena, CA 91125 

(Dated: July 2, 2008) 

We consider a class of inflationary models with small oscillations imprinted on an otherwise 
smooth inflaton potential. These oscillations are manifest as oscillations in the power spectrum of 
primordial perturbations, which then give rise to oscillating departures from the standard cosmic 
microwave background power spectrum. We show that current data from the Wilkinson Microwave 
Anisotropy Probe constrain the amplitude of a sinusoidal variation in the inflaton potential to have 
an amplitude less than 3 x 10~^. We anticipate that the smallest detectable such oscillations in 
Planck will be roughly an order of magnitude smaller, with slight improvements possible with a 
post-Planck cosmic-variance limited experiment. 

PACS numbers: 98.80.-k 
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I. INTRODUCTION 

Cosmic microwave background (CMB) experiments 
continue to be consistent with the simplest predictions 
of inflationary models, even as the data become increas- 
ingly precise 1] . Constraints on the amplitude and spec- 
tral index of primordial density perturbations, and on the 
amplitude of the inflationary gravitational-wave back- 
ground, can now be used to constrain the parameter 
space of the inflaton potential. In most of the current 
analyses the inflaton potential is parameterized in terms 
of its amplitude V{(j)) and first and second derivatives, 
V'[(j)) and y"(0). 

Given the recent advances in data quality, as well as the 
improvements anticipated with forthcoming experiments 
(e.g., the Planck satellite [2|, to be launched within the 
next year), it is worth asking whether the data can be 
used to study more complicated forms of the potential. 
For example, Adams et al. Q] showed that supergravity- 
inspired models may give rise to inflaton potentials with 
a large number of steps. Each step corresponds to a 
symmetry-breaking phase transition in a field coupled to 
the inflaton. The inflaton mass then changes suddenly 
when each transition occurs. These steps are responsi- 
ble for unusual inflaton dynamics, often represented as 
a hybrid inflation model 4], and they will create oscil- 
lating features in the primordial power spectrum. Os- 
cillations can also be directly imprinted on the infla- 
tonpotential itself, due to some trans-Planckian physics 
[3, [y, 0| • Finally, some other mechanisms may create fea- 
tures in the primordial power spectrum [8i]. There may 
also be empirical motivations to consider more compli- 
cated potentials, as several CMB analyses suggest that 
the CMB power spectrum may be better fit by primor- 
dial power spectra with features than by smooth power 
spectra H, [13, [HI, [3 ■ 

In this paper, we consider a class of inflationary models 
that feature periodic oscillations imprinted on a smooth 
inflaton potential, which then give rise to oscillations in 
the primordial power spectrum. We look for these os- 



cillations in the Wilkinson Microwave Anisotropy Probe 
(WMAP) data, and then determine the smallest oscil- 
lation amplitude that will be probed with forthcoming 
experiments. Our work is somewhat similar to that in 
Ref. [13], which considers oscillations in the primordial 
power spectrum that arise from a step in the inflaton 
potential, but differs in that we consider wiggles in the 
inflaton potential itself. There is also related work in 
Ref. [1J| , which considers oscillations in the CMB power 
spectrum from a rapid phase transition during inflation. 
However, the work most closely related to ours is that in 
Ref. [1^, which considers oscillations in the inflaton po- 
tential in natural-inflation models. While they focused 
on constraints from existing data, we forecast also the 
detectability with future measurements. 

The plan of this paper is as follows. Section [ll] in- 
troduces the model and discusses the calculation of the 
power spectrum. Section IIIII presents numerical results 
for the primordial, matter, and CMB power spectra. Sec- 
tion IIVI presents results of our search for oscillations in 
the WMAP data. Section |V] discusses the forecasts for 
the detectability of oscillations in future experiments; 
we consider here both Planck and a post-Planck cosmic- 
variance limited experiment. We summarize and provide 
some concluding remarks in Section IVTl 



II. THE OSCILLATING MODEL 

In order to study the effects of small oscillations in 
the inflaton potential, we begin with a simple base infla- 
tionary model that is consistent with current data. We 
consider the simplest such model, namely a quadratic 
inflaton potential, Vq{4>) = ^m^cp'^, with m ~ ^/Stt x 
10~^ Mpi. This potential has a corresponding CMB am- 
plitude A,, = 1.2 X 10"^ for a ACDM model. Once nor- 
malized with the best-fit WMAP5 data, we get Ug = 0.96 
for a pivot scale k^, = 0.002 Mpc~^, consistent with cur- 
rent constraints [Iq . 

We then superimpose on this smooth potential a sinu- 



soidal fluctuation to give a potential of the form 



F(0) = ^^m'cf? 
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parameterized by an amplitude of oscillations and a pa- 
rameter (3 that characterizes the frequency. The ampli- 
tude a is assumed to be small so that the inflaton does 
not get stuck in one of the local minima introduced in 
the potential by the oscillations; we discuss the precise 
constraint later. We choose the phase (5 = and ex- 
plain below why our results will be the same for different 
values. 

The homogeneous dynamics are dictated by the Fried- 
mann equation for the scale factor and the inflaton equa- 
tion of motion. The Friedmann equation (in units where 
SttG = c — h = 1) for a universe containing a scalar field 
(j){t) with potential V{<j>) is 



3H^ = 3 



= l^' + v{^), 



and the scalar-field equation of motion is 
4> + 3HtP = -^. 



(2) 



(3) 



We solve these coupled differential equations for the scale 
factor a{t) and scalar field (j){t) numerically. The solution 
is insensitive to our choice of initial conditions for 0(t), 
as the solution exhibits an attractor behaviour if the field 
begins high enough in the potential |18l |. 

We then turn our attention to the perturbations. We 
express the power spectrum V-jz of the primordial curva- 
ture perturbation with the horizon-crossing approxima- 
tion [i3,[Ii, 



Vn{k) = [1 - 2(2C + l)eH + 2Cm] 



2-k\q 



(4) 



where C = -2 + ln2 + fe ~ -0.73, with h the Euler- 
Mascheroni constant. The right-hand side is evaluated at 
k = aH, and the Hamilton-Jacobi slow-roll parameters 
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The next step is to relate a value of (/> to a comoving 
wavenumber k that crosses the horizon at that value of 
(p. To do so, we note that the number of e-foldings of 
inflation between a time t and the end of inflation is 



N{t) = In 



a (tend) 



and in terms of fc, it is 



N{k) ~ 55 - In 
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(6) 



(7) 



where oq is the scale factor today (which we choose to be 
flo = 1), and Ho = (/i/3000) Mpc"^ with h ~ 0.72. Note 
that the uncertainty around the 55 e-foldings of inflation 
at horizon crossing is about 5, and comes mainly from 
the uncertainty in the reheating process [l9|. 

The distance scale A that exits the horizon varies 
roughly as the exponential of the change Scf) in the infla- 
ton (j). Hence an oscillation in in the inflaton potential 
should give rise to oscillations in the logarithm of the 
wavenumber k in the primordial power spectrum, and 
thus in the logarithm of the CMB multipole moment £. 

We then calculate the temperature and polarization 
(TT, EE, and TE) spectra for the model using the CAMB 
code [201, using the current best-flt parameters^ for a 
ACDM model from the WMAP 5-year results [l6|. We 
neglect the B-mode polarization, as it comes about either 
through gravitational waves or gravitational lensing and 
is always small compared with the E-mode polarization. 



III. POWER SPECTRA 

We begin by showing in Fig. [T] the primordial power 
spectrum, as well as the present matter power spectrum, 
for the standard smooth potential and for an oscillat- 
ing potential with parameters chosen to be [a,/3] = [5 x 
10"'', 3 X 10~^]. The oscillations in the power spectra are 
nearly sinusoidal in log k with almost constant ampli- 
tude. The corresponding CMB angular power spectrum, 
shown in Fig. [2l reflects this behaviour. The mapping 
from the three-dimensional matter power spectrum to the 
two-dimensional CMB power spectrum slightly smooths 
the wiggles. In this example, the frequency (3 of the in- 
flationary potential's oscillations has been chosen so that 
one of the primordial oscillation peaks lines up with the 
first acoustic peak. However, this coincidence will not be 
generic. 



IV. A SEARCH IN WMAP DATA 

Given the CMB predictions discussed above, it is 
straightforward to search the existing CMB data from 
WMAP for these oscillations, and we have carried out a 
Markov Chain Monte Carlo (MCMC) analysis to do so. 
Fig. [3] presents the results of this analysis, assuming that 
all other cosmological parameters are known. For each 
point, we run four Markov chains. The results provide an 
indication of the best-fit amplitude a as a function of the 
assumed frequency (3. These amplitudes are comparable 
to their standard errors, and so we conclude that there is 
no evidence for oscillations in the WMAP data. Instead, 
we infer an upper limit a < 3 x 10~^ from WMAP. We 
note that a corresponding MCMC analysis has already 
been performed in Ref. [l3[ for the WMAP3 data for a 
similar potential. They find a precision on the amplitude 
of oscillations of the order of 10~^ for WMAP3, consis- 
tent with our results. 
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FIG. 3: The best-fit ampUtude a and its standard error, con- 
sidering WMAP5 data, for our range of frequencies /3. 



V. DETECTABILITY 

We now estimate the smallest oscillation amplitude a 
that will be detectable with future experiments. To do 
so, we first consider for simplicity the temperature power 
spectrum only. We suppose that each multipole moment 
£ can be measured with a standard error 



k [ h/Mpc ] 
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FIG. 1: Primordial power spectrum (top) and matter power 
spectrum (bottom) of the smooth inflaton potential (solid) 
and oscillating potential (dashed). The oscillating-potential 
parameters are [a,/3] = [5 x 10~*,3 x 10"'^]. The amplitude is 
chosen to be large to clearly show the effect of the oscillations. 



where 
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{Nfet) 
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is the contribution from the detector noise, AT is the 
detector noise per angular-resolution element, and Qu is 
the beam width. We then estimate the error on a by 21 1 
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FIG. 2: The CMB power spectrum corresponding to the mod- 
els shown in Fig.[T] The WMAP5 data are superimposed [la |. 
The error bars include both the cosmic variance and instru- 
mental noise. 
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and we choose an amplitude a — 5x 10^^ to compute the 
Ci derivatives. This estimate assumes that all other cos- 
mological parameters are known, and as such, provides 
an optimistic estimate. However, the true value, obtained 
by marginalizing over all other cosmological parameters, 
will probably not be too much worse, as there are no 
strong degeneracies between these oscillations and any 
other cosmological parameters. The acoustic oscillations 
oscillate in £, while these oscillate in log£, and so they 
should not be strongly degenerate. 

In order to improve our results, we then include the 
polarization and temperature-polarization power spectra 
as well. In that case, the generalization of the expression 
for the smallest detectable oscillation takes the form [221 
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FIG. 4: From the top to the bottom, by pairs, the smallest de- 
tectable amplitude cja as a function of /3, for WMAP, Planck 
and a cosmic-variance limited experiment. The upper part of 
each pair includes the temperature data only, and the lower 
one the polarization data as well. Note that the temperature 
and temperature-polarization curves for WMAP and Planck 
are effectively degenerate and thus appear to be one curve. 



for A = TT, EE, TE where [^-'^]aa' are elements of the 
inverse of ^, the covariance matrix; its elements are given 
inRef. dJ. 

In Figs. [3] and [3 we plot the smallest detectable am- 
plitude a a, from these analyses, as a function of /3, for 
WMAP (in Fig. 2] only), Planck, and a hypothetical 
cosmic- variance limited experiment. For WMAP, we sim- 
ply consider the measured uncertainties from the 5-year 
results [ia]. For Planck, the forecast is done considering 
the three most sensitive temperature channels, of spec- 
ifications similar to the HFI channels of frequency 100 
GHz, 143 GHz, and 217 GHz 0. The intensity sensi- 
tivities (detector noise) of these channels are taken as 
6.8 jjlK, 6.0 /iK, and 13.1 /iK respectively, correspond- 
ing to the values quoted for two complete sky surveys. 
These are average sensitivities per pixel, where a pixel is 
a square whose side is the FWHM extent of the beam. 
The FWHMs (beam widths) of these channels are given 
as 9.5 arcmin, 7.1 arcmin, and 5.0 arcmin, respectively. 
The composite noise spectrum for the three temperature 
channels is obtained by inverse-variance weighting the 
noise of individual channels [23, |2J| ■ For polarization we 
take only one channel, the 143 GHz channel, of sensitivity 
11.5 ^K, and FWHM 7.1 arcmin. Finally, for the cosmic- 
variance limited experiment, the fractional sky covered is 
taken to be 0.8 for all £, and we use simulated data out 
to an ^max of 2000, as for Planck. 

The statistical analysis presented in this Section is in- 
tended to forecast the uncertainty on parameters mea- 
sured by future experiments. It is nevertheless re- 
assuring that this forecast, when applied to WMAP, 
recovers roughly the limits obtained from the detailed 
MGMG analysis. Fig. |4] illustrates the advantage of 
Planck over WMAP for an oscillation search. While the 
smallest detectable amplitude is 0(10"^) for WMAP, it 
is O(10~^) for Planck. Looking at Planck more closely. 



FIG. 5: Magnification of Fig. |4] for Planck and the cosmic- 
variance limited experiment. 



in Fig. [5l we find a smooth variation of Ua except at 
some frequencies, such a.s (3 — 0.03 and /? — 0.04 for 
instance. These deviations arise from correlations be- 
tween the acoustic peaks and the primordial peaks, com- 
ing from the inflaton potential, in the CMB power spec- 
trum. When varying the frequency, the bumps thus cre- 
ated are aligned differently with the acoustic oscillations. 
For example, the choice of /3 = 0.03 considered in Fig. [2] 
is particular as many bumps are well aligned with the 
acoustic peaks. This effect is however much less impor- 
tant at high frequency, as our results suggest. 

The temperature-only and temperature-polarization 
curves coincide for both WMAP and Planck. This fol- 
lows because the polarization amplitude is much smaller 
and neither WMAP nor Planck will measure polarization 
to the cosmic- variance limit. The contribution to the to- 
tal signal-to-noise from polarization in these experiments 
is thus small. The degeneracy is broken in the cosmic- 
variance limited experiment. In this case, the precisions 
with which the temperature and polarization power spec- 
tra can be measured are roughly the same, resulting in 
roughly a v2 improvement to the oscillation sensitivity. 

Figs, m and [5] indicate that a higher frequency of os- 
cillations (i.e. a smaller /3) allows a smaller detectable 
oscillation amplitude. This can be understood by con- 
sidering Eq. pH)) . whose result is mainly governed by the 
difference between the chaotic and the oscillating curve, 
in terms of amplitudes. For a given amplitude of oscilla- 
tions in the CMB power spectrum, but different frequen- 
cies, the difference between these two curves is almost 
the same (it would be exactly the same if the slope of 
the chaotic spectrum was constant). However, in our 
analysis, we are considering a given amplitude a in the 
inflaton potential itself. When the frequency of oscilla- 
tions in the inflaton potential increases, so does that in 
the CMB power spectrum, but the corresponding ampli- 
tude in the CMB spectrum also increases (rather than 
staying constant). As such, for a given a, the difference 
between the chaotic and oscillating curve increases when 
(3 decreases, and hence CTq decreases. 



We can now also justify our neglect of the phase S in- 
troduced in Section [ni Its effect would simply be to shift 
the deviations on the smooth curves. In Section |lll we 
assumed that a was sufficiently small that inflation was 
not interrupted, and we can now justify that assumption. 
The constraint on a from WMAP is already a < 3x10^^. 
In the case of the most critical values considered in 
our analysis, being [a,/3] = [5 x 10~^,5 x 10"'^], we get 
jgmax^ |^max|] ^ ^^ ^ 10"^ 2 X 10"^ and [|?yg''''|, Iv'S"'']] ^ 
[2 X 10^^, 2 X 10^''], well within the slow-roll regime. 



VI. CONCLUSIONS 

The purpose of our paper has been to study whether 
features in the inflaton potential, beyond the standard 
slow-roll parameterization, can be detected with cur- 
rent and forthcoming CMB experiments. To do so, we 
have considered the possibility that the inflaton potential 
varies sinusoidally with the inflaton. A sinusoidal varia- 
tion of V{(l)) with (j) induces a variation in the matter and 
CMB power spectra that is sinusoidal in log A; and log^, 
respectively. 

Our analysis of the WMAP data indicates that the 
upper bound on a, the oscillation amplitude, is roughly 
3 X 10~^. To derive this bound, we argued that the de- 
generacy between a and other cosmological parameters 
is weak, and then assumed in the data analysis that the 
other parameters were fixed. We anticipate that a more 



complete analysis, including marginalization over other 
parameters, will weaken this bound, but only slightly. 
We then showed that the sensitivity of Planck to the os- 
cillation amplitude will be greater, relative to WMAP, 
by roughly an order of magnitude. For both Planck and 
WMAP, the constraint comes primarily from the tem- 
perature. Planck might be improved upon slightly by a 
cosmic-variance limited experiment. The advantage of 
such an experiment is that the additional information in 
the polarization may then be fully capitalized upon. 

On the theory side, we used the horizon-crossing ap- 
proximation to calculate the power spectrum. The work 
of Ref. [1^1 suggests that numerical integration of the 
perturbation equations will result in a slight suppression 
of the oscillation amplitude, relative to that we have ob- 
tained, but certainly by no more than a factor of two. 
If so, then our bounds may be accordingly weakened. 
We leave investigation of these data-analysis and theory 
questions for future research. 
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